Mesoscopic Fluctuations of the Loschmidt Echo 
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We investigate the time-dependent variance of the fidelity with which an initial narrow wavepacket 
is reconstructed after its dynamics is time-reversed with a perturbed Hamiltonian. In the semiclas- 
sical regime of perturbation, we show that the variance first rises algebraically up to a critical time 
t c , after which it decays. To leading order in the effective Planck's constant h e ff, this decay is given 
by the sum of a classical term ~ exp[— 2Af], a quantum term ~ 2h e s exp[— Ft] and a mixed term 
~ 2exp[— (F + \)t]. Compared to the behavior of the average fidelity, this allows for the extraction 
of the classical Lyapunov exponent A in a larger parameter range. Our results are confirmed by 
numerical simulations. 

PACS numbers: 74.40. +k, 05.45.Mt, 03.65.Yz 



Fluctuations of a physical quantity often contain more 
information than its average. For example, quantum sig- 
natures of classical chaos are absent of the average den- 
sity of states, but strongly affect spectral fluctuations Q. 
In the search for such signatures, another approach has 
been to investigate the sensitivity to an external pertur- 
bation that is exhibited by the quantum dynamics Q. 
Going back to Ref. 0, the central quantity in this ap- 
proach is the Loschmidt Echo 0, the fidelity 



M{t) = \(^ \exp^Ht}c^p[-iH Q t}\^ )\' 



(1) 



with which an initial quantum state i/'o is reconstructed 
after the dynamics is time-reversed using a perturbed 
Hamiltonian, H = Hq + eV (we set K = 1). This ap- 
proach proved very fruitful, however, most investigations 
of M(t) (which we will briefly summarize below) consid- 
ered the properties of the average fidelity M(t), either 
over different ipo, or different elements of an ensemble of 
unperturbed Hamiltonians H (having for instance the 
same classical Lyapunov exponent A) and/or perturba- 
tion V. Curiously enough, the variance a 2 (M) of the 
fidelity has been largely neglected so far. The purpose 
of this article is to fill this gap. We will see th at the 
variance er 2 (M) has a much richer behavior than M(t), 
allowing for the extraction of A in a larger parameter 
range, and exhibiting a nonmonotonous behav ior with a 
non-self-averaging maximal value a(t c )/M(t c ) ~ 1. 

We first su mmarize what is known about the average 
fidelity M(t) in quantum chaotic systems. Three regimes 
of perturbation strength are differentiated by three en- 
ergy scales the energy bandwidth B of Hq, the golden 

rule spreading F = 2Tre 2 \(4>a \V\4>^)\ 2 /A of an eigen- 



state 



/To) 



of Hq over the eigenbasis {</> a } of H, and 
the level spacing A = Bh c g (h e g = v d /£l is the effec- 
tive Planck's constant, given by the ratio of the wave- 
length volume to the system's volume). These three 
regimes are (i) the weak perturbation regime T < A, 
with a typical Gaussian decay M(t) ~ exp(— S 2 < 2 ), 



\(i> \V 2 \iP ) - (V'ol^o) 2 ), S 2 * TAh 



cir 



HQ 

(corrections to this Gaussian decay have been discussed 
in Ref. Q), (ii) the semiclassical golden rule regime 



A < r < B, where the decay i s expo nential with a rate 
set by the smallest of T and A, M(t) ~ exp[— min(r, X)t] 
Q10 01 and (iii) the strong perturbation regime T > B 
with another Gaussian decay M{t) ~ exp(— B 2 t 2 ) 
This classification is based on the scheme of Ref, 
which relates the behavior of M{t) to the local spec- 
tral density of eigenstates of Hq over the eigenbasis of 
H 0, Accordingly, regime (ii) corresponds to the 
range of validity of Fermi's golden rule, where the lo- 
cal spectral density has a Lorentzian shape 0, 0, Hcj . 
Quantum disordered systems with diffractive impurities, 
on the other hand, have been predicted to exhibit golden 
rule decay oc exp[— Tt] and Lyapunov decay cx exp[— At] 
in different time intervals for a single set of parameters 
[l^ . It is also worth mentioning that regular systems 
exhibit a very d ifferen t behavior, where in the semiclas- 
sical regime (ii), M(t) decays as a power-law 0] (see also 
Ref. [lj])- Finally, while in chaotic systems the averaging 
procedure has been found to be ergodic, i.e. considering 
different states Vo is equivalent to considering different 
realizations of H or V, the Lyapunov decay exists only 
for specific choices where ipo has a well defined classical 
meaning, like a coherent or a position state |3. ITlllTRIlfl |. 

Investigations beyond this qualitative picture have 
focused on crossover regions between the regimes (i) 
and (ii) @ and deviations from the behavior (ii) ~ 
cxp[— min(r, X)t] due to action correlations in weakly 
chaotic systems [13. Ref. [3 provides the only analyti- 
cal investigation of fluctuations of M(t) to date. It shows 
that, for classically large perturbations, r ^> B, M(t) is 
dominated by very few exceptional events, s o that a typ- 
ical ^o's fidelity is better described by exp[ln(M)], and 
that M(t) does not fluctuate after the Ehrenfest time 
tE = A _1 | ln[ft c ff]|. We will see that these conclusions do 
not apply to the regime (ii) of present interest. While 
some numerical data for the distribution of M(t) in the 
weak perturbation regime (i) were presented in Ref. [Tflj l . 
we focus here on chaotic systems and investigate the be- 
havior of a 2 (M) in the semiclassical regime (ii). 

We first follow a semiclassical approach along the lines 
of Ref. p|. We consider an initial Gaussian wavepacket 
M^o) = (™ 2 )- d/4 exp[ipo • (r' - r ) - \r' - r \ 2 /2v 2 }, 
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and approximate its time-evolution by 



r2 



(r|expH#oi)|Vo>= / dr' J2 A'f°(r, r' ; i)^o(r ), 



a 



1/2 



exp[iSf°(r,r ;t)-MT/i s /2]. 



The semiclassical propagator is expressed as a sum over 
classical trajectories (labelled s) connecting r and r' in 
the time t. For each s, the partial propagator contains 
the action integral (r, r ;i) along s, a Maslov index 
/is, and the determinant C s of the stability matrix [2lj . 
We recall that this approach allows to calculate the time 
evolution of smooth, localized wavepackets up to alge- 
braically long times oc Q(h~£) > t E (with a > 0)0. 
The fidelity then reads, 



M(t) 



dn / dr' / dr'i Mr'oWoW) 



(3) 



£jr£°(r lj r' ;t) [<( ri ,r' ';i)] a 



We want to calculate M 2 (t). Squaring Eq. <EJ, we see 
that M 2 {t) is given by eight sums over classical paths and 
twelve spatial integrations. Noting that ipo is a narrow 
Gaussian wavepacket, we first linearize all eight action 
integrals around ro, 

S s (r, r ;f) ~ S s (r, r ; t) - (r - r ) ■ p.. (4) 

We can then perform the Gaussian integrations over the 
eight initial positions r , r ' and so forth. In this way 
M 2 {t) is expressed as a sum over eight trajectories con- 
necting ro to four independent final points Tj over which 
one integrates, 

-.4 8 

M 2 (t) = / Y[drj cxp[i($ ffo - <£ H - ttM/2)] 



3 = 1 s;;i=l 



d/4 



I^ 1 / 2 - expH^/2) ,(5) 



where we introduced M = X^ =0 (— ^Y(t J, s2i+i ~ Ms 2 i+ 2 ) 
and 6p Si = p Si - p . 

The expression of Eq. Q is schematically described 
in Fig. Classical trajectories are represented by a full 
line if they correspond to Hq and a dashed line for H, 
with an arrow indicating the direction of propagation. In 
the semiclassical limit S s 3> 1 (we recall that actions are 
expressed in units of h), Eq. © is dominated by terms 
which satisfy a stationary phase condition, i.e. where the 
variation of the difference of the two action phases 

= S?°(r 1 ,v ;t)-S? 3 °(v ,r 2 ;t) 

+ ^ (r4,ro;i)-^ (r ,r 3 ;i), (6a) 

= S^(r ,r i; i) -S£(r 2 ,ro;f) 

+ S£(r 0l r 4 ;t) -fi£(r 3 ,r ;t), (6b) 
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Figure 1: Diagrammatic representation of the squared fidelity 
M 2 (t). 



has to be minimized. These stationary phase terms are 
easily identified from the diagrammatic representation as 
those where two classical trajectories s and s' of oppo- 
site direction of propagation are contracted, i.e. s = s' , 
up to a quantum resolution given by the wavelength v 
|23l |. This is represented in Fig. [51 by bringing two lines 
together in parallel. Contracting either two dashed or 
two full lines allows for an almost exact cancellation of 
the actions, hence an almost perturbation-independent 
contribution, up to a contribution arising from the finite 
resolution v with which the two paths overlap. How- 
ever when a full line is contracted with a dashed line, 
the resulting contribution still depends on the action 
SS S = —e J s V(q(t),t) accumulated by the perturbation 
along the classical path s, spatially parametrized as q(t). 

Since we are interested in the variance a 2 (M) = M 2 —M 
(this is indicated by brackets in Fig. |2j we must subtract 

——2 

the terms contained in M corresponding to independent 
contractions in each of the two subsets (si, S2, S3, S4) and 
(S5, S6, sj, ss). Consequently, all contributions to a 2 (M) 
require pairing of spatial coordinates, |r.; — rj| < u, for 
at least one pair of indices i,j = 1,2,3,4. 

With these considerations, the four dominant contri- 
butions to a 2 (M) are depicted on the right-hand side of 
Fig. The first one corresponds to si = S2 — S7 = s$ 
and S3 — Si ~ s,5 = se, which requires ri ~ r3, r2 ~ r4. 
This gives a contribution 



2</ 



J d ri dr 3 J2c 2 sl 



(7) 



x exp[-2^5pf 1 + i5® sl ]Q(v - |ri - r 3 | 



where S1> Sl = e f* dt'VV[q(t')][q Sl (f) - q S7 (t')\ arises 
from the linearization of V on s = s\_2 — s' = S7^ $.\vi\. 
and q Sl (?) lies on s\ with q(0) = r and q(t) = ri. In 
Eq. the integrations are restricted by |r x — r 3 | < v 
because of the finite resolution with which two paths 
can be equated (this is also enforced by the presence 



3 




a 2 {M) — a 2 e' 2Xt + 2ae- xt e~ Tt + 2h cS e- Tt Q(t - t E ) -f h 2 eS Q(t - t E ) 

Figure 2: Diagrammatic representation of the averaged fidelity variance a 2 (M) and the three time-dependent contributions 
that dominate semiclassically, together with the contribution giving the long-time saturation of a 2 (M). 



of i5$ s as we will see momentarily). For long enough 
times, t ^> t* , the phases 5$> s fluctuate randomly and 
exhibit no correlation between different trajectories [2(i| . 
One thus applies the Central Limit Theorem (CLT) 



(exp[i<5$ s 



exp[-(«S$2)/2] ~ exp[-e 2 /di(VV(0) 



VV(i)) |ri — r3| 2 /2A]. After performing a change of inte- 
gration variable J dr J2 S C-s — J dp and using the asymp- 
totic expression C s ~ (m/t) d exp[— At] [2J, one gets 



(T 2 = a 2 exp[— 2 At], 



2 t 2 /dr(Vy(0) • VV(t)) 



d/2 



(8a) 
(8b) 



The second dominant term is obtained from s\ = s 2 — 
s 7 = Sg, S3 = Si and S5 = sg, with ri ~ r3, or equiv- 
alently Si = S2, S7 = ss and S3 = S4 ~ S5 = sg with 
r 2 — r 4 . Therefore this term comes with a multiplicity 
of two, and one obtains 

a\ = 2(^) 2 ^/^r 3 £C 2 1 

x exp[-2^ 2 ( 5p 2 1 + t5^ Sl \e(iy - |n - r 3 |)\ 

x (/ dr2 C °* c M-v 2 Sp 2 S3 + MS B9 ]\ , (9) 

again with the restriction |ri — i"3| < i^. To calculate the 
first bracket on the right-hand side of Eq. QJ, we first 
average the complex exponential, assuming again that 
enough time has elapsed so that actions are randomized. 
The CLT gives (exp[i55 S3 ]) = exp(-±{5S% a )) with 



«> = 



dt dF(y[q(i)]y[q(?)]). (10) 
Jo 



Here q(t) lies on S3 with q(0) = r and q(t) = r 2 . In 
hyperbolic systems, correlators typically decay exponen- 
tially fast, 



(V[q(t)]V[q(t')])^exp[-r]\t-t'\], 



(11) 



with an upper bound on r\ set by the smallest positive 
Lyapunov exponent [24| • One thus obtains (<5S' 2 3 ) = Tt. 
Usually T cx e 2 is identified with the golden rule spread- 
ing of eigenstates of H over those of H 0, 0] ■ It is dom- 
inated by the short-time behavior of (V r [q(t)]V r [q(0)]). 



We stress however that for long enough times, (SSg ) oc t 
still holds to leading order even with a power-law decay 
of the correlator (V[q(t)]V[q(?)]) oc \t-t'\-^, provided 
77 is sufficiently large, 77 > 1. We note that similar ex- 
pressions as Eq. (II fl relating the decay of M to time 
integrations over the perturbation correlator have been 
derived in Refs. 0, El using a different approach than 
the semiclassical method of Ref. used here. Further 
using the sum rule 



(^ 2 A) C 



dr 



C s exp[-^<5p 2 ] =1, 



one finally obtains 

cr| = 2aexp[-At] exp[— Tt]. 



(12) 



(13) 



The third and last dominant time-dependent term 
arises from either si = S7, s 2 = sg, S3 = S4, S5 = sq 
and ri ~ r 3 , or s\ — s 2 , s 3 = s 5 , s 4 = s 6 , S7 = s 8 and 
r 2 — r 4 . It thus also has a multiplicity of two and reads 

V —\ N dv 1 dv 2 dv z dv i Y J c s 1 Cs 2 C S3 C S5 
x exp[-v 2 (5p 2 Sl + <5p 2 2 + <5p 2 3 + 5p 2 J] 



x 



exp[i(6S S3 - 5S S5 )] e(i/-|n-r 3 |)>. (14) 



The integrations, again, have to be performed with |ri — 
r 3 | < v. We incorporate this restriction in the calculation 
by making the ergodicity assumption, setting 



dr 1 dr 2 dr 3 dr 4 . ..®(v - |n - r 3 |)) 



H eS ( / dr 1 dr 2 dr 3 dr 4: . . .}<d(t - t E ), 



(15) 



which is valid for times larger than the Ehrenfest time 
[111 (for shorter times, t < tE, the third diagram on the 
right-hand side of Fig. goes into the second one). One 
then averages the phases using the CLT to get 



of = 2h cS exp[-Tt\e{t-t E )- 



(16) 



Subdominant terms are obtained by higher-order con- 
tractions (e.g. setting r 2 ~ r 4 in the second and third 
graphs on the right hand-side of FigEJ. They either de- 
cay faster, or are of higher order in H e s, or both. We only 



4 



discuss the term which gives the long-time saturation at 
the ergodic value a 2 {M) ~ ft 2 ff . For t > tE, there is 
a phase-free (and hence time-independent) contribution 
with four different paths, resulting from the contraction 
si = s r , s 2 = s 8 , s 3 = 85, s 4 = s 6 , and ri ~ r 3 , r 2 ~ r 4 . 
Its contribution is sketched as the fourth diagram on the 
right-hand side of Fig. It gives 

(2 \ '2d „ 
V -\ [J dcidxa^2C n C aa (17) 

x expl-^^p^ + Sp 2 J}e(u - |n - r 3 |)) 2 . 

From the sum rule of Eq. (112L and again invoking 
the long-time ergodicity of the semiclassical dynamics, 
Eq. Qjj}, one obtains the long-time saturation of <r 2 (M), 



0.5- 

p(M) r 

03- 



0.1- 




-4-2 0_ 2 4 
(M - M)/CT 

Figure 3: Distribution P(M) of the fidelity computed for 10 4 
different ^0 for N = 32768, <5A" = 5.75 ■ 10" 5 (i.e. T « 0.09), 
at times f = 25, 50, 75 and 100 kicks. 



al = hl s Q{t-t E ). 



(18) c 2 (M) has a nonmonotonous behavior, i.e. it first rises 



Note that for t < tE, this contribution does not exist by 
itself and is included in o\, Eq. JBJ. 

According to our semiclassical approach, the fidelity 
has a variance given to leading order by the sum of the 
four terms of Eqs. ©, (ES), OS) and OBI 

cr s 2 c = a 2 exp[-2Ai] + 2a exp[-(A + T)i] (19) 

+2fi cfr cxp[-r£]e(i - t E ) + h 2 cS e(t - t E ). 

Eq. illHfl is the central result of this paper. We see that 
for short enough times, i.e. before ergodicity and the 
saturation of M{t) ~ H e g and cr 2 (M) ~ h 2 s is reached, 
the first term on the right-hand side of iflfljl will dominate 
as long as A < r. For A > Y on the other hand, <j 2 {M) 
exhibits a behavior cx exp[— (A + T)t] for t < tE, turning 
into (x h e g ex.p[—rt] for t > tE- Thus, contrary to M, 
tr 2 (M) allows to extract the Lyapunov exponent from the 
second term on the right-hand side of Eq. IjliJH even when 
A > r. Also one sees that, unlike the strong perturbation 
regime r » 5 ^3], M(t) continues to fluctuate above 
the residual variance ~ fi, 2 ff up to a time ~ r~ 1 |ln/i G ff| 
in the semiclassical regime B > T > A. For T <C A, 
lnftcffl 3> t E and M{t) fluctuates beyond tE- 
The above semiclassical approach breaks down at short 
times for which not enough phase is accumulated to mo- 
tivate a stationary phase approximation j23|. To get 
the short-time behavior of er 2 (M), we instead Taylor 
expand the time-evolution exponentials exp[±ii2>o)t] = 
i 2 /2 + ... + 0(Hf m t 5 ). The resulting 



1 ± iH (0) t 



expression for <r 2 (M) contains matrix elements such as 
(■0o|-ff( a o - ) |V'o), a = 1,2,3,4, which one then calculates 
using a Random Matrix Theory (RMT) approach l2fl l 
for the chaotic quantized Hamiltonian iTm) 0, @> llflj . 
Keeping non-vanishing terms of lowest order in t, one 

has a quartic onset o 2 {M) ~ {TA - £2 2 )t 4 for t < E _1 , 
with E a = [e 2 ((^ \V 2 \iP ) - {MV\^?)] a/2 - RMT gives 
(E 4 — E 2 ) cx (TB) 2 , with a system-dependent prefactor 
of order one. From this and Eq. I|19J1 one concludes that 



at short times, until it decays after a time t c which one 
can evaluate by solving cr 2 c (t c ) = (r_B) 2 t 4 . In the regime 
B > r > A one gets 



\TBJ 



1 A (—} 1/2+d 1 
\TBJ 2 + d 



(20) 



and thus 



a 2 (t c ) ~ (TByi^) 



2 (ao\ i/2+d 

\tb) 



4A 



2 + d \TBJ 



(21) 



We explicitely took the t-dependence a(t) = aoi _d into 
account. We estimate that ao cx (T\)~ d / 2 (obtained by 
setting the Lyapunov time equal to few times the time 
of flight through a correlation length of the perturbation 
potential, as is the case for billiards or maps), to get 
a 2 {t c ) cx {B/X) 2d l 2+d > 1. Because < M{t) < 1, this 

value is however bounded by M (t c ). Since in the other 
regime T < A, one has er 2 (^ c ) ~ 2h ofi [l-(2^ off ) 1 / 4 1 /r/B] 
we predict that <r 2 (t c ) grows during the crossover from 
T <C A to T > A, until it saturates at a non-self- averaging 
value, a(t c )/ M(t c ) sa 1, independently on % e s and B, 
with possibly a weak dependence on V and A. 

We conclude this analytical section by mentioning that 
applying the RMT approach to longer times reproduces 
Eq. O with A — > 00 plj. This reflects the fact that 
RMT is strictly recovered for = only. 

To illustrate our results, we present some numerical 
data. We based our simulations on the kicked rotator 
model with Hamiltonian |23l 



Hn 



2 



K cos x 8{t — n) . 



(22) 
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We concentrate on the regime K > 7, for which the 
dynamics is fully chaotic with a Lyapunov exponent 
A = hi[AT/2]. We quantize this Hamiltonian on a torus, 
which requires to consider discrete values pi = 2irl/N and 
xi = 2nl/N, I = 1, ...TV, hence h cS = 1/N. The fidelity 
(QJ is computed for discrete times t = n, as 



M(n) = \(M(UZK) n (Uo) n \i>o)f 



(23) 



using the unitary Floquet operators Uq = 
cxp[— ip 2 /2h c fi] exp[— iKq cosx/h c s] and Usk hav- 
ing a perturbed Hamiltonian H with K = Ko + SK. 
The quantization procedure results in a matrix form 
of the Floquet operators, whose matrix elements in 
x— representation are given by 



(Uo) 



1 



/./' 



N 



cxp[z- 



.ir(l 



V\2 



N 



exp[ 



2wl' 
-— — cos — — 

2tt N 1 



The local spectral density of eigenstates of Usk over those 
of Uo has a Lorentzian shape with a width T cx (5K/h e s) 2 
(there is a weak dependence of T in -Ko) in the range 
5 = 27r>r>A = 2n/N). This is illustrated in the 
inset to Fig. 

Numerically, the time-evolution of tpo m the fidelity, 
Eq. J 2 .'ill , is calculated by recursive calls to a fast-Fourier 
transform routine. Thanks to this algorithm, the matrix- 
vector multiplication Uq^k^o requires O(TVlnTV) oper- 
ations instead of 0(N 2 ), and thus allows to deal with 
very large system sizes. Our data to be presented below 
correspond to system sizes of up to N < 262144 = 2 18 
which still allowed to collect enough statistics for the cal- 
culation of a 2 (M). 

We now present our numerical results. Fig. shows 
the distribution P(M) of M(t) in the regime F < A for 
different times. It is seen that even though P(M) is not 
normally distributed, it is still well characterized by its 
variance. A calculation of o 2 (M) is thus meaningful. 

We next focus on <r 2 in the golden rule regime with 
r«A, Data are shown in Fig. ||| One sees that a 2 (M) 
first rises up to a time t c , after which it decays. The max- 
imal value cr 2 (t c ) in that regime increases with increasing 
perturbation, i.e. increasing T. Beyond t c , the decay of 
a 2 is very well captured by Eq. (II fill . once enough time 
has elapsed. This is due to the increase of a 2 (t c ) above 
the self- averaging value cx h e g as T increases. Once the 
influence of the peak disappears, the decay of a 2 (M) is 
very well captured by erf given in Eq. I|lf>|) . without any 
adjustable free parameter. Finally, at large times, a 2 (M) 
saturates at the value given in Eo.l|18jl. 

As SK increases, so does V and a 2 (M) decays faster 
and faster to its saturation value until T > A. Once 
T starts to exceed A, the decay saturates at exp(— 2Xt). 
This is shown in Fig. which corroborates the Lya- 
punov decay of a 2 (M) predicted by Eqs. ©. Note 
that in Fig. the decay exponent differs from the Lya- 
punov exponent A = ln[if/2] due to the fact that the 
fidelity averages (C s ) cx (exp[— At}) ^ exp[— (X)t] over 
finite-time fluctuations of the Lyapunov exponent [T3 |. 




Figure 4: Variance <r 2 (Af) of the fidelity vs. t for weak 
r < A, N = 16384 and 10 s ■ 8K = 5.9, 8.9 and 14.7 (thick 
solid lines), N = 4096 and SK = 2.4 ■ 1CT 4 (dashed line) 
and N = 65536 and SK = 1.48 ■ 1(T 5 (dotted-dashed line). 
All data have Ko — 9.95. The thin solid lines indicate the 
decays = 2h cB exp[-rt], with V = 0.024(5^ ■ iV) 2 (there is no 
adjustable free parameter). The variance has been calculated 
from 10 3 different initial states ipo- 




Figure 5: Variance cr 2 (Af) of the fidelity vs. t in the golden 
rule regime with T > A for TV = 65536, Ko = 9.95 and 
SK £ [3.9 ■ 10" 5 , 1.1 • 1(T 3 ] (open symbols), and N = 262144, 
Ko = 9.95, SK = 5.9 ■ 1(T 5 (full triangles). The solid line is 
cx exp[— 2Ait], with an exponent Ai = 1.1, smaller than the 
Lyapunov exponent A = 1.6, because the fidelity averages 
(cxp[— Xt]} (see text). The two dashed lines give ft 2 ff = TV -2 . 
In all cases, the variance has been calculated from 10 3 
different initial states tpo- 



At long times, a 2 {M) saturates at the ergodic value 
a 2 (M,t — > oo) = h 2 s , as predicted. Finally, it is seen 
in both Figs. 0] and El that t c decreases as the perturba- 
tion is cranked up. Moreover, there is no A^-dependence 
of <7 2 {t c ) at fixed T. These two facts are at least in qual- 
itative, if not quantitative, agreement with Eq. 1(201) . 

The behavior of a 2 {t c ) as a function of T is finally 
shown in Fig. El First we show in the inset the behavior 
of the local spectral density 



E 



(24) 



of eigenstates {</>i ^} (with quasienergy eigenvalues e a ) 
of Uq over the eigenstates {<fi a } (with quasienergy eigen- 




io 3 w 2 10 

r/B 



Figure 6: Maximal variance cr 2 (t c ) as a function of F/B, 
for K = 10.45, N = 4096, 16384, 65536 and 262144 (empty 
symbols) and K = 50.45, N = 16384 (full circles). The 
variance has been calculated from 10 3 different initial states 
4>o- Inset: local spectral density of states p(e) of eigenstates 
of an unperturbed kicked rotator with Ko = 12.56 over the 
eigenstates of a perturbed kicked rotator with K = K + 8K, 
SK = 5 ■ 10~ 3 . System sizes are N = 250 (diamonds), 
N = 500 (circles) and N = 1000 (squares). The solid lines 
are Lorentzian with widths F « 0.0125, 0.05 and 0.0124 in 
agreement with the formula L = 0.024 (SK ■ N) 2 . 



values e a ) of Usk- As mentioned above, p(e) has a 
Lorentzian shape with a width given by T « 0.024((5J4T • 
A^) 2 . Having extracted the N— and 5K— dependence of 
r, we next plot in the main part of Fig. U3 the maximum 
a 2 (t c ) of the fidelity variance as a function of the rescaled 
width T/B of p(e). As anticipated, cr 2 (t c ) first increases 
with T until it saturates at a value > 0.1, independently 
on hog, T or A, once r^B. These data confirm Eq. i'lll 
and the accompanying reasoning. Note that once T ex- 
ceeds the bandwidth B, p(e) is no longer Lorentzian, and 
the decay of both M(t) and a 2 (M) is no longer exponen- 
tial 0|. 

In conclusion we have applied both a semiclassical and 
a RMT approach to calculate the variance a 2 (M) of the 
fidelity M(t) of Eq. QJ. We found that a 2 (M) exhibits a 
nonmonotonous behavior with time, first increasing alge- 
braically, before decaying exponentially at larger times. 
The maximum value of a 2 (M) is characterized by a non- 
self-averaging behavior when the perturbation becomes 
sizable against the system's Lyapunov exponent. 
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